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Research  Objectives 

Particle  and  droplet  size  distributions,  being  parameters  of 
fundamental  importance,  should  be  priority  measurement  objectives 
for  intelligent  sensors  in  next  generation  propulsion  systems. 
Unfortunately  there  are  a  number  of  problematic  scientific  issues 
impeding  the  development  of  laser  light  scattering  particle 
sizing  instruments  capable  of  on-line,  autonomous,  and  self- 
diagnosing  operation  in  hostile  environments.  The  objective  of 
this  research  effort  is  to  contribute  to  the  scientific  knowledge 
base  necessary  to  characterize  and  hopefully  extend  the 
applicability  of  laser  diffraction  instruments  under  these 
adverse  conditions.  Both  the  novel  specific  approaches  discussed 
below  and  the  emphasis  on  concepts  relevant  to  intelligent 
instruments  and  sensors  make  this  research  project  unique. 
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1.  Inverse  Scattering  Algorithms 
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Particle  sizing  using  laser  (Fraunhofer)  diffraction  is  an 
inverse  scattering  problem.  The  governing  equation  for  single 
scattering  by  large  particles  is  a  Fredholm  integral  equation: 
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where:  a  =  vd/X  is  the  particle  size  parameter;  n(a)  is  the 
unknown  particlp  size  distribution;  and  1(0  the  scattering 
intensity  which  can  be  measured  at  some  number  of  scattering 
angles.  Unfortunately  the  presence  of  n(a)  inside  the  integral 
of  Eq.  (1)  makes  the  problem  ill-posed  in  the  sense  that  small 
variations  in  the  measurements  of  1(0  (e.g.  due  to  noise)  can 
cause  large  changes  in  the  resulting  n(a).  The  maximum 
information  obtainable  about  n(a)  is  dependent  on  the  noise 
levels  in  the  measurements,  the  sampling  or  detection  scheme,  the 
stability  or  robustness  of  the  inversion  method,  the  actual  size 
distribution,  and  the  allowable  level  of  uncertainty.  Ideally, 
then,  a  laser  diffraction  system  would  have  variable  detector 
geometry,  inversion  scheme,  and  size  resolution  which  could  all 
be  adjusted  at  run  time  depending  on  the  conditions  encountered. 


The  objective  of  this  part  of  the  research  is  to  identify 
and  impact  the  underlying  scientific  issues  which  presently 
impede  the  needed  quantum  improvements  in  efficiency  and 
robustness  of  inversion  schemes  for  intelligent  laser  diffraction 
sensors.  Dynamic  adaptive  control  of  the  detector  configuration, 
synergism  of  the  inversion  scheme  and  detector  geometry  are  major 
focus  areas. 

2.  Multiple  Scattering. 

A  major  challenge  for  future  research  on  in-situ  optical 
techniques  for  particle  sizing  is  to  extend  application  of  these 
methods  into  the  regime  of  optically  thick  media.  The 
theoretical  framework  for  particle  size  measurements  under 
multiple  scattering  conditions  is  rather  tenuous,  even  for  the 
direct  problem  of  calculating  the  radiation  transfer  properties 
of  an  arbitrary  particulate  medium.  Even  less  developed  are 
analytical  methods  for  the  inverse  problem,  i.e.  obtaining 
particle  size  distribution  information  from  measured  optical 
scattering  properties  of  media  where  multiple  scattering  is 
significant . 


The  objective  of  this  part  of  the  research  is  to  understand 
the  process  of  multiple  scattering  as  it  might  occur  in  future 
applications  of  intelligent  particle  size  distribution  sensors 
and  develop  inversion  schemes  which  can  operate  in  such  an 
environment.  A  minimum  requirement  is  that  the  presence  of 
multiple  scattering  be  diagnosed  by  the  instrument  to  ensure  that 
erroneous  data  not  be  used  in  control  algorithms,  a  situation 
which  could  result  in  a  catastrophic  failure.  Our  objective  is 
to  surpass  that  and  develop  efficient  and  robust  algorithms  which 
can  actually  extract  useful  particle  size  information  from 
Fraunhofer  diffraction  measurements  in  multiple  scattering 
environments . 
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3.  Laser  Beam  Deflection/Steering 


Another  problem  which  has  already  limited  the  usefulness  of 
laser  diffraction  particle  sizing  instruments  in  propulsion 
systems  research  is  that  of  beam  deflection  by  refractive  index 
gradients  in  the  optical  path.  These  gradients  are  due  to 
spatially  nonuniform  temperature  and/or  species  concentrations 
and  at  relatively  large  scales  have  the  same  effect  as  an  array 
of  randomly  dispersed  and  time-varying  lenses  which  cause  the 
probe  beam  to  move  around  at  the  detection  plane.  This  has  been  a 
severe  problem  since  information  on  the  largest  droplets  in 
typical  sprays  is  concentrated  at  small  scattering  angles  which 
are  of  the  same  order  as  the  beam  deflection  angles  encountered. 

The  objective  for  this  research  topic  is  to  investigate 
concepts  which  could  enable  a  laser  diffraction  instrument  to 
operate  under  conditions  where  beam  deflection  is  significant.  A 
minimal  requirement  is  for  on-line  detection  of  beam  steering  so 
that  an  intelligent  instrument  would  at  least  know  that  the  data 
taken  under  these  conditions  is  suspect.  A  more  important 
objective  is  to  develop  strategies  to  permit  a  diffraction  system 
to  autonomously  perform  real-time  correction  for  beam  steering. 

Status  of  the  Research  Effort 

Significant  progress  has  been  made  in  each  of  these  areas 
over  the  duration  of  this  grant  and  we  foresee  no  significant 
deviations  from  the  original  approach  proposed  for  the  research. 
The  sections  below  describe  the  technical  achievements  made  in 
these  areas  during  the  indicated  12  month  time  period. 

1.  Inverse  Scattering  Algorithms 

The  inverse  scattering  problem  for  laser  diffraction 
particle  sizing  instruments  has  similarities  with  all  systems 
requiring  a  deconvolution  of  experimental  data.  The  inverse 
scattering  problem  (i.e.  determination  of  n(a)  in  Eq.  (1)  from 
measured  1(0))  can  be  solved  by  either  (i)  direct  integral 
transform  methods  or  (ii)  indirect  numerical  quadrature 
approaches.  In  the  past  year  we  have  concentrated  on  the 
generally  overlooked  integral  transform  techniques  which,  by 
virtue  of  the  direct  nature,  have  unique  potential  for  very  fast 
solutions  as  would  be  required  by  an  on-line,  intelligent  sensor. 
The  general  integral  transform  solution  to  Eq.  (1)  is: 

n(«)  =  c/«z  J  hx(«0)  X(0)  d 0  (2) 

>,,n 

where:  c  is  a  constant,  hx  is  the  transform  kernel,  X  is  some 
function  of  the  experimental  data,  and  the  specific  forms  of  X 
and  hx  are  coupled.  The  most  straightforward  solution  to  Eq.  (2) 
obtains  an  X(0)  proportional  to  the  derivative  of  1(0). 
Integration  of  Eq.  (2)  by  parts  produces  other  functions  X  (and 
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associated  functions  hx  )  which  involve  integrals  of  1(6).  We  have 
initiated  a  thorough  study  of  the  integral  transform  formulations 
applicable  to  the  inverse  Fraunhofer  diffraction  problem  and 
identified  four  unique  methods.  The  methods  are  identified  by  the 
primary  author  as:  (i)  Chin  (1955),  (ii)  Shi f rin-Kolmakov  I 
(1967),  (iii)  Shi f r i n-Kolmak ov  IT  (1967),  and  (iv)  Petrov 
(1976/1985).  Only  for  the  method  of  Chin  had  complete  studies  of 
inversion  performance  been  previously  made. 


In  order  to  establish  a  framework  to  evaluate  existing  and 
new  transform  methods  we  have  developed  a  set  of  benchmark 
problems.  These  benchmark  data  sets  permit  numerical 
s t ud i es /s imul a t i ons  of  the  performance  of  various  techniques 
under  controlled  and  realistic  conditions.  In  general,  there  are 
six  major  parameters  that  will  affect  the  inversion  of  laser 
diffraction  data.  They  are  (a)  minimum  scattering  angle  6m  .  ,  (b) 
maximum  scattering  angle  0nax,  (c)  angular  resolution  &6 ,  (d) 
noise  level,  (e)  out  of  range  particles,  and  (f)  numerical 
techniques  for  integration  and  differentiation.  As  a  part  of 
this  research  we  have  studied  inversion  results  obtained  for 
assumed  bi-modal  and  tri-modal  log  normal  particle  size 
distributions  as  a  function  of  these  six  parameters. 
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An  important  property  of  an  inversion  algorithm  is  immunity 
to  experimental  measurement  error  or  noise.  Figure  2  indicates 
the  results  of  a  Chin  inversion  of  synthetic  scattering  data 
perturbed  by  several  forms  of  synthetic  noise.  The  results  are 
reasonably  good,  even  for  the  rather  large  levels  of  noise  used 
in  the  simulations.  The  other  integral  transform  methods  perform 
in  a  manner  similar  to  that  shown  in  Fig.  2  and  the  results  are 
documented  in  the  publications  listed  below. 

2.  Multiple  Scattering 

A  solution  to  the  inverse  multiple  scattering  could  be 
obtained  if  an  instrument  could  ascertain  which  photons  reaching 
the  detector  plane  had  been  scattered  once  and  only  once.  This 
single  scattering  data  could  then  be  inverted  using  methods  based 
on  Eq.  (1)  as  discussed  above.  Unfortunately  we  have  been  unable 
to  suggest  a  method  for  performing  this  photon  segregation 
directly.  An  indirect  method  we  are  developing  involves  viewing 
the  effect  of  the  multiple  scattering  medium  as  a  matrix 
operation  such  that: 

u  I  t  =  M  •  S6  ingle  ( 3 ) 

where:  Sa u ] t  and  SB 4 n  j e  are  vectors  with  each  element 
representing  the  scattering  signal  collected  for  one  of  the  ring 
detectors  of  Fig.  1  and  the  subscripts  indicate  multiple  and 
single  scattering  respectively;  and  M  is  a  square  matrix  which 
represents  the  multiple  scattering  characteristics  of  the  medium. 
S1(llt  is  the  scattering  signature  which  would  be  measured 
directly  by  a  laser  diffraction  system.  Mj ^  represents  the 
efficiency  by  which  the  medium  of  interest  redistributes  that 


single  scattered  energy  which  would  have  reached  the  jth  detector 
(under  optically  thin  conditions)  into  the  ith  detector.  If  the 
matrix  M  could  be  measured,  then  the  effective  single  scattering 
signature  Sg . n  j e  could  be  determined  by  inverting  M: 

Single  =  •  Sau:t  (4) 

The  size  distribution  could  then  he  determined  using  a  classical 
inversion  of  the  effective  single  scattering  signature  using 
standard  inverse  Fraunhofer  algorithms. 

We  have  developed  and  are  studying  a  method  to  measure  the 
matrix  M  on-line.  Consider  a  diffraction  reticle  of  Fig.  3  which 
could  generate  a  hollow  cone  of  scattered  light.  This  synthetic 
scattering  signature  would  represent  an  input  vector  Slin  )(  with 
only  one  non-zero  element.  After  this  passed  through  the  medium 
the  diffraction  signature  at  each  detector  due  solely  to  the 
input  scattering  cone  would  be  proportional  to  a  corresponding 
element  in  a  single  column  of  M  as  seen  in  Eq.  (3).  Figure  3  is  a 
computer-generated  hologram  designed  to  produce  the  hollow 
scattering  cone.  A  series  of  scaled  holograms  like  Fig.  3,  one 
for  each  detector,  provides  a  series  of  scattering  vectors  S, 
each  with  a  different  non-zero  element.  In  that  way  M  can,  in 
theory,  be  characterized  on-line.  Our  experiments  to  investigate 
the  performance  of  the  diffraction  hologram  are  continuing. 

3.  Laser  Beam  Deflection/Steering 

The  probe  laser  beam  from  a  Fraunhofer  diffraction 
instrument  will,  in  general,  be  continuously  deflected  as  it 
traverses  a  medium  of  nonhomogeneous  refractive  index  as  would  be 
encountered  in  a  combustor.  This  causes  the  transmitted  beam  and 
the  associated  diffraction  pattern  which  is  centered  around  it 
to  be  directed  to  off-center  points  in  the  detector  plane  of  Fig. 
1.  The  very  high  intensity  in  the  transmitted  beam,  which 
typically  would  be  more  than  lOOx  greater  than  the  intensity 
scattered  into  the  small  angle  inner  rings  of  the  detector, 
severely  hampers  the  measurements  at  small  angles.  One  possible 
solution  we  are  investigating  is  to  measure  the  instantaneous 
position  of  the  probe  beam  at  the  detector  plane  and  dynamically 
configure  as  ring  detector  centered  about  that.  We  have  obtained 
a  position  detector  as  shown  in  Fig.  1  and  find  that  the  time- 
response  is  adequate.  Research  into  the  remaining  problem  of 
configuring  a  detector  in  real-time  is  continuing. 

Technical  Presentation 
1 .  E .  D .  Hirleman 

"Calibration  Techniques  for  Laser  Diffraction  and  Single 

Particle  Counting  Spray  Sizing  Instruments"" 

Fifth  ASTM  Symposium  on  Pesticide  Formulations  and 

Application  Systems,  Kansas  City,  MO.,  November,  1984. 
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Technical  Publications 


1.  E.  D.  Hirleman  and  L.  G.  Dodge 

"Laser  Diffraction  Measurements  of  Sprays  and  Reticles: 
Accuracy  and  Interpretation" 

Western  States  Section/Combustion  Institute  Spring  Meeting, 
San  Antonio,  TX,  April  1985,  Paper  85-6. 

2.  E.  D.  Hirleman,  E.  Hovenac,  and  R.  F.  Ide 
"Calibration  and  Sample  Volume:  Characterization  of  PMS 
Optical  Array  Probes" 

ICLASS  ’85  Proceedings,  Paper  85-84,  The  Institute  of 
Energy,  London,  England,  July,  1985. 

3.  E.  D.  Hirleman  and  L.  G.  Dodge 

"Performance  of  Laser  Diffraction  Drop  Size  Analyzers" 
ICLASS  *85  Proceedings,  Paper  85-64,  The  Institute  of 
Energy,  London,  England,  July,  1985. 


Personnel 

Prof.  E.  Dan  Hirleman  -  Associate  Professor  of  Mechanical  and 
Aerospace  Engineering.  Principal  Investigator 

Joseph  H.  Koo  -  Research  Assistant  and  Ph.D.  student.  Mr.  Koo’s 
dissertation  will  come  out  of  this  work  and  his  Ph.D.  should  be 
awarded  in  1986. 

New  Discoveries 

1.  Computer- generated  Hologram  for  Diffraction  Reference 

Signature 

One  new  development  resulting  from  this  research  is  the 
concept  of  a  computer-generated  hologram  to  produce  a  diffraction 
reference  signature  which  has  only  a  single  scattering  order 
(i.e.  scattering  at  a  single  angle).  A  first  generation  of  this 
concept  has  been  produced  and  a  photographic  reproduction  is 
shown  in  Fig.  3.  This  concept  is  an  important  advance  in  the 
future  development  of  inverse  Fraunhofer  diffraction  schemes 
which  may  work  in  multiple  scattering  environments.  The  concept 
also  has  important  implications  for  improvements  in  on-line 
calibration  of  laser  diffraction  instruments. 
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